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Abstract 



The magnetic field is found to inhibit the Ekman suction. The width of the Ekman layer is reduced 



^sO ■ In this paper we present axisymmetric nonlinear simulations of magnetized Ekman and Stew- 

artson layers in a magnetized Taylor-Couette flow with a centrifugally stable angular-momemtum 
profile and with a magnetic Reynolds number below the threshold of magnetorotational instability 

i 

with increased magnetic field normal to the end plate. A uniformly-rotating region forms near the 

1*3 

O ' outer cylinder. A strong magnetic field leads to a steady Stewartson layer emanating from the 

£>y junction between differentially rotating rings at the endcaps. The Stewartson layer becomes thin- 

Qh< ner with larger Reynolds number and penetrates deeper into the bulk flow with stronger magnetic 

| field and larger Reynolds number. However, at Reynolds number larger than a critical value ~ 600, 

> 

£T) ' axisymmetric, and perhaps also nonaxisymmetric, instabilities occur and result in a less prominent 

^J- ' 

■ Stewartson layer that extends less far from the boundary. 
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I. INTRODUCTION 



The history of Taylor- Couette flow dates back to the 19th century. To measure viscosity, 
Couette n studied Sows between rotattng concentric cylinders. Rayleigh's staMity criterion 
was introduced in 1916 during his study of cyclones. Taylor [2J extended it by including 
viscosity, and made quantitative predictions of instability in Couette flow. If the cylinders 
were infinitely long, the steady-state laminar solution would be the ideal Taylor-Couette 
state: 

Qo{r) = a+±, (1) 

in which a = (f^?"! — ^i r i)/( r i — r i) an d b = rfr|(fii — ^2) /(^I — r i)> where T\ and r 2 
are the radius of the inner and outer cylinder, and Qi and Q 2 are the angular velocity of 
the inner and outer cylinder respectively. Rayleigh's stability criterion states that in the 
unmagnetized and inviscid limit, such a flow is linearly axisymmetrically stable if and only 
if the specific angular momentum increases outwards: ab > 0. 

The study of magnetized Taylor-Couette flow began much later. Velikhov and Chan- 
drasekhar 4j discovered that a vertical magnetic field may destabilize the flow, provided that 
the angular velocity decreases outward, Vl\ < which today is called magnetorotational 
instability (MRI). In ideal magnetohydrodynamics (MHD), the instability takes place with 
an arbitrarily weak field {5], Bl. Experiments on magnetized Couette flow aiming to observe 
MRI have been performed [7|, [8|, but MRI has never been conclusively demonstrated in 
the laboratory. Some other experiments have been proposed or are still under construc- 



tion 



EH, 



11 



12|. The experimental geometry planned by most groups is a magnetized 



Taylor-Couette flow: an incompressible liquid metal with density p, kinematic viscosity v 
and magnetic resistivity rj confined between concentric rotating cylinders, with an imposed 
axial and/or toroidal background magnetic field sustained by currents external to the fluid. 

The challenge for experimentation is that liquid-metal flows are very far from ideal on 
laboratory scales. While the fluid Reynolds number Re = Qiri(r 2 — ri)/u can be large, the 
corresponding magnetic Reynolds number Re m = fiiri(r2— Ti)/r] is modest or small, because 
the magnetic Prandtl number Pr m = u/rj ~ 10~ 5 — 10 -6 in liquid metals. Standard MRI 
modes will not grow unless both the rotation period and the Alfven crossing time are shorter 
than the timescale for magnetic diffusion. This requires both Re m > 1 and S > 1, where S = 
Va{t 2 — T\)/r] is the Lundquist number, in which Va = B z0 /^/fi p is the Alfven speed and B z0 



2 



is the imposed axial magnetic field. Therefore, Re > 10 6 and fields of several kilogauss must 
typically be achieved. Hollerbach and collaborators have discovered that MRI-like modes 
may grow at much reduced Re m and S in the presence of a helical background field, a current- 

nn n 

free combination of axial and toroidal field [131 114J|. Though Stefani et al. [8J have claimed to 
observe this helical MRI (HMRI) experimentally, we explained the experimentally measured 
wave patterns to be transiently amplified disturbances launched by viscous boundary layers 



rather than globally unstable modes 



161 ] . We also questioned the relevance of this helical 



MRI to astrophysics by showing that this new mode is stable for a Keplerian rotation profile 
by WKB analysis in a narrow-gap geometry (see Sec. II. A of Liu et al. and by linear 
calculations in a wide-gap geometry (see Sec. II. B of Liu et al. [l^]). Recently Riidiger 

n 

and Hollerbach [17|, Priede et al. 
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have reported that this new mode is unstable in the 
inductionless limit with for some boundary conditions. Under the parameters used in the 
Riidiger and Hollerbach , Priede et al. 18] (Rm m = S = Pr m = 0, but with finite Re and 
Hartmann number Ha = Va[ti ~ r i)l \fw)-> the authors are indeed taking the diffusivity to 
infinity rj — > oo. Note however that the combination Ha 2 /2Re = (Va) 2 / [ZQi]) , which is the 
Elsasser number A [20], is also finite. The authors consider Ha and Re to be constant as 
Pr m — > 0; thus if we think of Q and v as fixed, then the Alfven speed must scale like y/rj as 
i] — > oo. So, the authors are considering a limit in which the Alfven speed is infinitely larger 
than the rotation speed but poorly coupled to the flow, whereas Liu et al. [ijj] were thinking 
of the resistive limit as one in which the field strength and rotation speed were held fixed as 
the diffusivity became infinite. The former limit is unlikely to be important in astrophysics. 
However, it might be achieved in a low-plasma-beta but highly resistive (weakly ionized) 
plasma. 

In view of the large Reynolds number, the Taylor-Proudman theorem suggests that the 
end plates should dominate the entire flow unless a very long cylinder were used (h/(r2 — 
ri) ^> 10 3 , where h is the height of the cylinders.) ]19]. The no-slip boundary condition 
on the end plates causes an imbalance between centrifugal and pressure forces and drives 
Ekman circulation. If the endcaps rotate rigidly with the outer cylinder, this circulation 
takes the form of inward flows along the endcaps, which turn vertically along the inner 
cylinder, converge at the midplane, and depart the cylinder in a radial jet j^]. This Ekman 
circulation, and especially the jet, transport angular momentum efficiently and reduce the 
free energy avai.aiffe for shear-driven instabilities Q. Both effect s are unfavorable for 
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laboratory demonstration of MRI. The Princeton MRI experimental apparatus has been 
constructed to minimize the circulation by the use of independently controlled split endcaps 



20 



21 



221 ] . Nevertheless the jump of the rotation speed at the junction of the rings extends 



some distance into the bulk Stewartson layer 



24j . however, the modification of the 



Stewartson layer by the axial magnetic field has to be studied. 



There has been research done on the MHD Ekman 
as they are sometimes called in the literatures) 
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avers 



has been little work aside from Pariev 



23, 



or Ekman-Hartmann layers 



28 



30 



31] . However there 



32] concerning the effect of finite differential rota- 



tion on Ekman 
(Hollerbach 
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ayers, or on magnetized Stewartson layers in cylindrical geometry 23j, [24 ] 



34J ] did excellent work on magnetized Stewartson layers, though in spherical 



geometry) . The latter issues remain poorly understood but play a big role in MRI experi- 



ments 



16l . |35| and have potential importance in geophysics and fluid dynamics. Hollerbach 



and Fournier [19] discussed the purely hydrodynamic (unmagnetized) steady results with 
the assumption of infinitesimal differential rotation, or a very tiny Rossby number, while 

our paper discusses time-dependent solutions with finite differential rotation. Szklarski and 

I 

Riidiger |35] presented results with finite differential rotation but without rings (similar to 
Sec. Ill of our paper), thus no Stewartson layer is present. That paper is a good contribution 
to Potsdam Rossendorf Magnetic Instability Experiment (PROMISE) and the discussion of 
Taylor-Dean flow is very insightful. This paper is one of the first to study magnetized 
Stewartson layers in cylindrical geometry. Understanding the role of the boundary layers, 
especially magnetized ones, is critical to the success of MRI experiments. 

It is known that Ekman circulation is significantly modified when the Elsasser number 



29] exceeds unity: 



A = B 2 J(8np V n) > 1 , (2) 

where Q is the characteristic rotation frequency, which we take equal to ^2- For gallium, 
A « 2.5(-B/Tesla) 2 (1000 rpm/O). Adopting the parameters used in the Princeton MRI 
experiment (TableJT]), Q = Q2 = 533 rpm, B = 5000 Gauss, immediately leads to A ~ 1.2. 
In the PROMISE experiment [8], A ~ 2.4. Hence magnetic modifications to the Ekman 
layer should be significant in both experiments. 



Here we report nonlinear simulations with the ZEUS-2D code 36j, [37], which is a time- 
explicit, compressible, astrophysical ideal magnetohydro dynamics (MHD) two-dimensional 
code, to which we have added viscosity and resistivity (with subcycling to reduce the cost 



of the induction equation) for axisymmetric flows in cylindrical coordinates (r,ip,z) 38]. 
The simulation domain mimics the Princeton MRI experiment (see Table |T] and Fig. [1]) ex- 



cept where stated explicit 
in Sec. II. D of Liu et al. 



y. The code adopts the magnetic boundary conditions introduced 



15l | (not the commonly used vertically pseudo-vacuum boundary 



conditions). All real flows are actually compressible; in an ideal gas of fixed total volume, 
density changes generally scale ~ M 2 when Mach number M = Vfi ow /V sonri d < 1. Incom- 
pressibility is an idealization in the limit M — > 0. An isothermal equation of state has 
been used with a sound speed chosen so that the maximum of M < 1/4. The techniques 
used here have been benchmarked analytically and compared with other codes in Liu et al. 



15l . |38| and verified experimentally in Liu et al. [161 ] . Burin et al. [21|. Note that in the 



simulations the magnetic diffusivity rj is fixed to the experimental value rj ~ 2, 000 cm 2 s 1 
(TableJT]), however the kinematic viscosity is varied for the purpose of extrapolation. Also 
as demonstrated by Goodman and Ji Q, the viscosity of liquid metals is so small as to be 
almost irrelevant to MRI, at least in the linear regime. 

In this present work, since we concentrate on magnetic Ekman and Stewartson layers, the 
rotation speed profile is chosen so that the system is MRI stable. We have found the MRI 
linear growth rates (7 s" 1 ) as a function of magnetic Reynolds number Re m and Lundquist 
number S (Fig. [2]) from a WKB analysis [U| with the same dimensions as in Table [J and 
fj, = f^/^i = 0.13325, which shows that for Re m < 10 the system should be MRI stable. 
All simulations presented in this paper are in this regime. 

This paper is organized as follows. We start in Sec. [II] by reviewing a one- dimensional ap- 
proximation to magnetized Ekman circulation above an infinite, uniformly rotating bound- 
ary. Two-dimensional effects are introduced in Sec. IHH but still with rigid endcaps. In 
Sec. IIV[ we divide each endcap into two independently rotating rings as in the Prince- 
ton MRI experiment, and study the dependence of the resulting Stewartson layer on the 
Reynolds and Elsasser numbers. Implications for the Princeton MRI experiment are dis- 
cussed in Sec. [V] 
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Figure 1: Computational domain for studies of magnetic Ekman layer. Region (I): Perfect conduct- 
ing inner cylinder, angular velocity infinitely long. (II): Liquid metal. (Ill): Perfectly insulating 
inner ring, Q3, extending to infinity; (IV): Perfectly insulating outer ring, O4, extending to infinity; 
(V): Perfectly conducting outer cylinder, infinitely long. Thin dash line: the midplane. B z is 
the initial background vertical uniform magnetic field. 



II. STANDARD MAGNETIC EKMAN LAYER WITH NEAR-UNIFORM ROTA- 
TION 



We begin with a problem considered by Gilman and Benton [25|]. The problem treated 
consists of an incompressible, viscous and resistive fluid above an infinite, flat and insulating 
boundary that rotates at angular velocity O = Qe z . Far from the boundary, the fluid rotates 
uniformly at Q — fi(l + e). A uniform magnetic field aligned with the rotation axis is 



imposed. In the amdysis of Gilman and Benton 25], an expansion in powers of e, together 
with von Karman similarity |39l . 1401 ]. leads to a solution that is exact to first order in e. In 
the limit that e < 1, increasing A results in a continuous transition between pure Ekman 
flow and a rotating analog of Hartmann flow. 

Here we sketch a modified steady state WKB analysis rather than an expansion in the 
von Karman similarity variables used by Gilman and Benton [25] . With the t and r depen- 
dence factored out, the linearized equations of motion reduce to inhomogeneous ordinary 
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Figure 2: MRI linear growth rates (7 s -1 ) as a function of magnetic Reynolds number Re m and 
Lundquist number S with dimensions as in Table HI n = 7.1 cm, r<i = 20.3 cm and h = 27.9 cm 
and fi = = 0.13325. The system is MRI stable if Re m < 10 regardless of S. 



differential equations with coefficients independent of z. Elementary homogeneous solutions 
of these equations exist with exponential dependence on z\ however, since there is an insu- 
lating endplate at z = 0, the wavenumber k n may be complex, and the final solution can be 
a linear combination of the elementary modes for different k n and one particular solution 
that matches the flow far from the boundary. The vertical magnetic boundary conditions 
require the fields to match onto a vacuum solution at the end plate. 
We seek a mode of the form 

8 

[v r ,v v , B r ,B v ] T = [0, Vqq, 0, 0] T + y^C n [ B r , n ,B^ n } T exp(ik n z). (3) 

n=l 

The first column vector on the right-hand side is the particular solution, which satisfies 
the boundary conditions at z = 00 but not at z — 0, where Voo is the velocity far away 
from the end plate. Each term in the sum above is the elementary solution corresponding 
to a particular wavenumber k n , with (v r n , . . . , B v n ) T a 4-component column vector; these 
elementary solutions are superposed with constant weights {C n }, which must be chosen to 
satisfy the boundary condition. The 8 values of the wavenumber {k n } are the roots of the 
steady-state dispersion relation: 

7 



Dimensions 


7*i = 7.1 cm 


r2 = 20.3 cm 


h = 27.9 cm 


Material Property 


p ~ 6.0 gem 3 


Tj ~ 2.0 x 10 3 cm 2 s 1 


Full Speed Run 


Q 1 /2ir = 4000 rpm 


0,2 /2ir = 533 rpm 


Q 3 /2tt = 1820 rpm 


^1^/217 = 650 rpm 


Rotation Profile used in Sec. IIIII 


n 1 /2ir = 500 rpm 


Q 2 /2ir = 66.625 rpm 


n 3 /2ir = 66.625 rpm 


Q 4 /2tt = 66.625 rpm 


Rotation Profile used in Sec. IIVI 


Oi/2vr = 500 rpm 


Q 2 /2ir = 66.625 rpm 


Q 3 /2ir = 227.5 rpm 


n 4 /2ir = 81.25 rpm 



Table I: Parameters used in the simulations 



k 4 [( V u) 2 k 4 + 2-qvVlk 2 + (Vi + 4^V)] = , (4) 

which follows from the linearized, homogeneous and axisymmetric Navier-Stokes and in- 
duction equations. Only the four nonzero roots of this equation are of interest since they 
determine the boundary-layer thickness. The eigenmodes corresponding to k = would 
modify the interior flow (z — > oo). The four nonzero roots are 

k * = Yl±™. (5) 

r/is v 

The two "acceptable" nonzero roots of Eq. [5l satisfying the boundary conditions, are k± = 
— (Jcr ± iki), where kn = <5 -1 as given by EqnEl so that 

v r = -V OQ e~ kRZ sin faz , (6) 
v 9 = Vcoil-e-*** cos k T z), (7) 
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where kj is related to &r by kf/k R = \/l + A 2 — A. Thus 

1 [ 1 — A/2 if A < 1 ; 

5 = S E &5 E x{ (8) 

V VA 2 + 1 + A 1/V2A ifA>l. 



Here 5 E = yv/Vt is the purely hydrodynamical Ekman-layer thickness with near-uniform 
rotation. Notably, the Elsasser number A (Eq. [2]) has nothing to do with v. Hence even if 
the boundary layer were turbulent, with an effective turbulent viscosity vj- and thickness 
increased by O^z/j/z/) 1 / 2 ], the magnetic field would be at least as consequential as in the 
laminar case. This assumes that turbulent magnetic diffusivity is negligible, as one might 
expect since the laminar value is large enough. One expects that A ^> 1 probably results in 
a more stable layer and pushes the onset of turbulence to larger Reynolds numbers. In the 
limit A — > oo, the thickness 5 — > ^yurj/VA. this is the Hartmann-layer thickness, which does 
not depend upon Q. 

The above theoretical results have been used to benchmark our code (Fig. [3]). The 
thickness of the Ekman layer 5 is the reciprocal of &r, which is deduced by fitting the 
simulated data at = (ri + r 2 )/2 using Eq. El The results agree well with the theoretical 
prediction (Eq. [8]). 

III. MAGNETIC EKMAN LAYER WITH DIFFERENTIAL ROTATION: (I) END- 
PLATES COROTATING WITH OUTER CYLINDER 

In contrast with the idealized case in Sec. [Til most Taylor- Couette experiments have 
a two-dimensional circulation driven by differences in the rotation of the inner and outer 
cylinders and the endcaps. In this section, we take the end plates to corotate with the outer 
cylinder, i.e., Q3 = = ^2- 

The Reynolds number based on the Ekman layer thickness is [32I ]: 

Res ^ r J^l^ Re y^3xl0 3 , (9) 

for full-speed runs of the Princeton MRI experiment (Table [TJ. The Ekman layer with 
uniform rotation as in Sec. [TTJ has two known instabilities, viscous and inflection point in- 
stabilities, both of which are axisymmetric. The viscous instability owes its existence to 



the perturbed Coriolis force 



411 ] while the inflection point instability is of the inviscid type. 
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Figure 3: The thickness of the Ekman layer 5 versus Elsasser Number A for Re = 1600, Re m = 5. 
fii/2vr = 1000 rpm, VL 2 /2k = 1000 rpm, n 3 /2ir = 1010 rpm, ^ 4 /2vr = 1010 rpm. n = 15 cm, 
r 2 = 35 cm and h = 20 cm. The data is measured at r = (r\ + r 2 )/2 = 20 cm. The dashed line 
is the theoretical result. The solid line is the one obtained from modified ZEUS-2D simulations. 
We have chosen larger n and r 2 than the ones of the Princeton MRI experiment to minimize the 
curvilinear effects and larger h IOSe to minimize the influence of the top endcap. 

n 

From Fig. 3 of Gilman [29] , the critical Reynolds number of the viscous and inflection point 
instabilities associated with this Ekman-Hartmann layer is in the range: 100 < Re$ c 1000 
for A « 1, at least for cases of near-uniform rotation as in Sec. [Ill Thus the boundary 
layer is turbulent for the full-speed runs of the Princeton MRI experiment. However in the 
simulations, the Reynolds number in the bulk is taken to be 6400, thus Res = 80, so that the 
boundary layer is laminar. Our discussion below is grounded on the equations of laminar 
flows. The magnetic Reynolds number in the boundary layer based on the thickness of the 
Ekman layer is defined as SjJ: 

Rms= 6 -^~^L, (10) 
V VRe 



where U is a characteristic speed. For Re = 6400 and Re m = 2.5 as in the simulations, 
Rms ~ 3.125 x 10~ 2 . Because Rms 1 and \u = (Q — f^/f^l ^ 1> the quantity 
1 + (l/2)r*duj/dr* = l/(2Qr)d(r 2 fl)/dr, which is the ratio of vorticity to rotation (r* is the 
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radius normalized by r 2 ), is ~ a/Q 2 > 0. The solution decays with oscillation asz-> oo, as 



for an unmagnetized Ekman layer. The modified Ekman layer thickness S is given by 



32|: 



S = <) L (n i r l fi, 1 (11) 

A2 + 1 + I r *f: + A 

Eq. [TT] reproduces Eq. [H]if there is no differential rotation (duj/dr* — > 0). Therefore the strong 
magnetic field causes the Ekman layer to become thinner even with differential rotation. It 
is worth emphasizing that the above derivation is based on a first order expansion in u <C 1. 

Our simulations with the parameters of Table [T] approach the regime n of the above linear 
theory except that: (1) the radial boundary condition is conducting rather than insulating 
(a magnetic Ekman layer with fully insulating boundaries on all sides is the next step for this 
problem and will be included in a forthcoming paper); (2) the flow profile far away from the 
end plate differs from the ideal Couette profile, though not greatly; (3) \(fl — f^A^I "C 1 
is not satisfied except near the outer cylinder. 

The simulations are analyzed at = (ri + r 2 )/2 to minimize two-dimensional effects due 
to the radial boundary. Since \(Q(rf) — f^)/^! = 1-08 is not small, some nonlinear effects 
neglected in the above linear analysis could be important. 

From Fig. HI we confirm that the axial magnetic field does reduce the Ekman layer 
thickness. The finite differential rotation cannot be neglected and modifies the linear Ekman 
layer, which is seen from the unmagnetized case, i.e. A = 0: The theoretical result predicts 
that the thickness of the Ekman layer with finite differential rotation is larger than the 
thickness of the Ekman layer with uniform rotation; in the numerical results, this is less 
clear. Though the simulated curve does not match the theoretical result very well, the 
agreement is as good as might be expected when a theory based on w < 1 is applied to 
simulations at u ~ 1. 

For Re = 6400, the final state is not steady even after at least five Ekman times te = 
hj '(z/k/2) 1 / 2 . Given a finite differential rotation, it is more appropriate to estimate Ekman 
time te via the epicyclic frequency, 



k = 2 ((rfrl - r*0?)/(r} - r») 



1/2 



which is the maximum frequency of small axisymmetric inertial oscillations inside the inviscid 
fluid. Typical (instantaneous) flow and field patterns are shown in Fig. [5j The poloidal flux 
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Figure 4: The thickness of the Ekman layer 8 versus Elssaser Number A for Re = 6400, Re m = 2.5. 
Parameters as in Table HI The data are measured at r = (r\ +r2)/2 = 13.7 cm. The dashed line is 
from the linear analysis (Eq. Hip . The solid line is obtained from modified ZEUS-2D simulations. 



and stream functions are denned so that 

V P = V r e r + V z e z = r~ x e^ x V $, B P = B r e r + B z e z = r _1 e^ x V (12) 

which imply V • Vp = and V • Bp = 0. Two Ekman cells are clearly visible. The 
flapping "jet" at the midplane due to the Ekman circulation breaks the vertical refelction 



symmetry of the system, resulting in a chaotic region around the midplane [42j . The poloidal 
flow circulation and toroidal field are small compared to the background toroidal flow and 
initial axial field respectively, 

max— ^— < 13% , max^- < 3% . 
riih B z0 

The most noticeable feature of the final state of the magnetic Ekman circulation is the 
presence of an area of solid body rotation near the outer cylinder (Fig. 0(d)) as in Holierbach 



33J . This area increases with the Elsasser number: the strong axial magnetic field squeezes 
the dynamically active area (Ekman cells) toward the inner cylinder (Fig. [6]). When A = 1.5, 
almost half of the liquid metal is rotating with the outer cylinder. This is due at least in part 
to the following two effects: (1) Larger axial magnetic fields suppress the Ekman circulation 

12 




Figure 5: Contour plots of final-state velocities and fields with uniformly rotating endcaps. Re = 
6400, Re m = 2.5 with B z0 = f500 Gauss (A = 1.09). Parameters as in Tabled (a) Poloidal flux 
function ^ (Gauss cm 2 ) (b) Poloidal stream function <E>(cm 2 s -1 ) (c) toroidal field (Gauss) (d) 
angular velocity f2 = r~ x V^ (rads -1 ). 

more thoroughly; (2) The axial Hartmann current turns towards the radial direction near 
the midplane, and couples with the axial magnetic field to produce an azimuthal Lorentz 
force, which tends to reduce the azimuthal velocity shear dVt/dz. Both effects reinforce 
Taylor-Proudman theorem near the inner cylinder, where there is a large velocity shear 



13 



between the bulk flow and the end plate. 
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Elsasser Number A 

Figure 6: AR/(r2 — r\) v.s. A. AR is the radial width of the poloidal circulating region. Re m = 2.5, 
Re = 6400 with endplates corotating with outer cylinder. Parameters as in TableUl From the radial 
profile of the azimuthal velocity iv,, AR is the radial gap between the inner cylinder and the place 
where the solid body rotation starts. Note that in the simulations the magnetic diffusivity rj is 
fixed to rj ~ 2,000 cm 2 s -1 , however the imposed axial magnetic field B z q is varied. 

IV. MAGNETIC EKMAN LAYER WITH DIFFERENTIAL ROTATION: (II) END 
PLATES SPLIT INTO TWO RINGS 

We have brought the computation closer to the experimental conditions by making the 
endcaps consist of two independently rotating rings as in Fig. [Hand TableUl The junction 
between these two rings lies at = (r x + r 2 )/2 = 13.7 cm. For Re = 6400, the final 
state is not steady. Typical (instantaneous) flow and field patterns are shown in Fig. [7J 
Two flapping "jets" due to the unsteady Stewartson layer, emanating from the junction of 
the rings at both endcaps, leads to a chaotic region localized there (Fig. [7J (b)), which is 
different from the case in Sec. IHIl in which the unsteady region is mainly near the midplane. 
The poloidal flow circulation and toroidal field are also small compared to the background 
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toroidal flow and initial axial field respectively, 

max^- < 4.3% , max|^ < 1.3%. 

These ratios are smaller than the ones discussed in Sec. IHIl which implies that the Ekman 
suction is reduced by splitting the endcaps into two differentially rotating rings. 




8 10 12 14 16 18 20 8 10 12 14 16 18 20 

R (cm) R (cm) 



Figure 7: Like Fig. El but with two differential rotating rings and A = 1.5. Parameters as in 
Table [H The Stewartson layer is located between the rings at = (n + ?"2)/2 = 13.7 cm and 
breaks the two big Ekman cells into eight cells. 
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The following observations can be made from Fig. [7J With rings the Stewartson layer 
is more apparent than Ekman circulation. The split endcaps break the two big Ekman 
cells found in Sec. II I II (Fig. [7J (b)). The four cells at intermediate radii are straightforward 
consequences of the Stewartson layer as discussed below. The direction of the circulation 
of the bottom four cells is opposite to the direction of the circulation of the corresponding 
upper cells. Hereafter we focus only on the upper half of the flow. 

The increase of the number of Ekman cells can be understood from Fig. [HJ The direction 
of the residual Ekman flow depends upon the angular velocity of the boundary relative 
to the interior, thus resulting in anti-clock-wise normal Ekman cells at r < 10.6 cm and 
13.7 cm < r < 18.2 cm and clock-wise abnormal Ekman cells elsewhere. 




8 10 12 14 16 18 20 
R [cm] 



Figure 8: Two independently rotating rings generate 8 cells, solid line, ideal Couette state; dashed 
line, rotation profile at the endcaps. Arrows indicate the radial flow directions near the endcaps. 

The magnetic field tends to reduce fluctuations in the final state at high Reynolds number. 
The Stewartson layer becomes more prominent with increasing A at fixed Re (Fig. [9]). This 
can be understood by considering the influence of the magnetic field on the stability of the 
Stewartson layer. The addition of an axial magnetic field (in the MRI stable regime) resists 
shear along the magnetic field lines and elongates the cells vertically so that they penetrate 
deeply into the fluid. More details are given in Sec. [V] 
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On the other hand, the Stewartson layer also becomes sharper as Re increases at fixed A 
From Fig. [10] (a), we infer the following scaling law: 

fi = 3.9 + 0.014Re - 57 . 
or 

This is somewhat consistent with the one-dimensional analyses of a purely hydrodynamic 
Stewartson layer, which show that a Stewartson layer consists of nested layers of outer thick- 



ness E 1 ^ and inner thickness £ 1/3 UM 
the idealizations used in the analyses 



where E = 1/ Re is the Ekman number. Considering 
241 ] and the complications in our two-dimensional sim- 
ulations, the agreement is as good as might be expected. 

We also observe that the Stewartson layer penetrates more deeply into the bulk flow with 
larger Re at fixed A, at least for Re < 400 and A = 1.5 in axisymmetry. This can be 
seen from Fig. [10] (b). The profiles deviate from the ideal Couette state more with larger 
Re [Re < 400). However at even higher Reynolds number (i.e., Re > 400), the Stewartson 
layer develops axisymmetric MRI/centrifugal instabilities at large axial wave numbers, which 
upon saturating result in a less prominent Stewartson layer, as Fig. [10] (b) (Re > 400). The 
layer is more localized near the endcaps in these cases. More details are given in Sec. [V] 

V. DISCUSSION 

Purely hydrodynamical (i.e., A = 0) experimental results show that the azimuthal ve- 
locity profi.e is auite smooth Q; no obvious Stewartson .aye r is observed at distances 
greater than ~ 1 cm from the bottom endcap [43|]. The effect of the velocity jump across 
the junction between rings is not as severe as in the simulations reported here and those 
by Hollerbach and Fournier [3]. This difference may be explained by various instabilities 
associated with the Ekman and Stewartson layers. 

As for the Ekman layer, the low Reynolds number used in the simulations leads to a 
laminar Ekman-Hartmann layer as discussed in Sec. IHIl The experimental boundary- layer 
Reynolds number Res = y/Re ~ 3 x 10 3 is larger than the critical value Re cr u ~ 10 3 given 



A ~ 1 [29(|, for the axisymmetric instabilities (viscous and inflection point instabilities) of 
Ekman-Hartmann layer with near-uniform rotation. And the vertical velocity shear due 
to the finite differential rotation in the experiment would result in the Kelvin-Helmholtz 
instability given a sufficiently high Reynolds number, which, however, could not be resolved 
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Figure 9: Azimuthal velocity cm s" 1 versus radius r at different heights with Re m = 2.5, 
Re = 6400, and the endcaps divided into two rings at = (n + r2)/2 = 13.7 cm. Parameters 
as in Table HI Solid line, ideal Couette state; long dashes, z = 1.33 cm (relative to the bottom 
endcap); dash dot, z = 2.79 cm; short dashes, z = 13.95 cm. (a): A = 0.38; (b): A = 1.5. 

by our axisymmetric simulations due to the same reason stated below. Therefore unstable 
Ekman layers are highly possible in the experiment. The layers may be smoothed by localized 
circulation and/or turbulence from these instabilities. This may account for the differences 
in the extent and prominence of the Stewartson layer between simulation and experiment. 

As for the Stewartson layer, at the junction of the rings, the outer ring rotates more slowly 
than the inner one (f^ < f2 3 ), hence d(r 2 Q 2 )/dr < across the junction. This radial shear 
could result in both the Kelvin-Helmholtz instability and Rayieigh centrifugal instability 
given a sufficiently high Reynolds number. Unfortunately our axisymmetric simulation could 



no 



tz instability since it is a toroidal nonaxisymmetric mode. 



44j , Hollerbach 45j , Hollerbach et al. 46( , Schaeffer and Cardin 



not resolve the Kelvin-Helm! 
However, see Friih and Read 

4?J for experimental and theoretical studies of such instabilities in other contexts. 

It is well known that surface tension at the interface between two fluids hinders the 
Kelvin-Helmholtz instability. In the real magnetized experiment, besides the instabilities 
discussed below, a less prominent Stewartson layer due to the Kelvin-Helmholtz instability 
would be resulted. In a homogeneous but magnetized fluid such as ours, magnetic field 
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Figure 10: (a) Time-averaged \d$l/dr\ of the Stewartson layer at z = 1.33 cm v.s. Reynolds number 
Re. (b) Time-averaged — ^o(^d)]/^o('"d)| with = (n + f2)/2 = 13.7 cm at the midplane. 
v.s. Re, where fio is the ideal Couette state defined in Eq. [TJ Re m = 2.5, A = 1.5, and the endcaps 
divided into two rings at = (n + ?"2)/2 = 13.7 cm. Parameters as in Table [Q Note that m 
the simulations the rotation speed profile is fixed to Q\/2-k = 500 rpm, f^/^vr = 66.625 rpm, 
^l^/2n = 227.5 rpm, 0,^/271 = 81.25 rpm, however the kinematic viscosity v is varied. 



tension may stabilize the Kelvin-Helmholtz instability 48|, |49|. However, the discussion 
of the nonaxisymmetric Kelvin-Helmholtz instability is beyond the scope of our current 
methods based on the axisymmetry and is not the main purpose of this paper. Similarly 
magnetic field tension may stabilize the Rayleigh's centrifugal instability 50( . We find that 
the short wavelength modes are stabilized before (i.e., at lower A) the long wavelength 
modes by performing a WKB stability analysis for the flows near the junction of the rings 
(Fig. [IT]) , where the Stewartson layer lies. Our analysis assumes axisymmetry, thus the 
Kelvin-Helmholtz instability is excluded. Following Ji et al. Ill], we have the dispersion 
relation: 



k 2 f)0 2 
[( 7 + uk 2 )^ + V k 2 ) + [k z V A f}% + k 2 ( 7 + r^ 2 ) 2 + £=-(k,V A 

hi omr 







(13) 



All variables have the same meanings as in Ji et al. Uj except: (1) the characteristic 
rotation speed fl is chosen to be v^3^4i (2) the dimensionless vorticity parameter, ( = 
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Figure 11: Growth rates (7 s _1 ) predicted by Eq. [13] as a function of Elsasser number A and 
vorticity parameter £ at the junction of rings, = (ri + r^jl with Re = 6400. Parameters as in 
Table HI Panel (a) the vertical mode number n = 1; panel (b) n = 2. 



(l/rn)d(r 2 n)/dr = 2 + 91nfi/<91nr is taken to be 2 + (r d /fi)(Afi/Ar), where Afi = fi 3 -fi 4 
and Ar is the radial thickness of the Stewartson layer; (3) the wave number k = \/k% + k%, 
where the axial wave number k z = nn/(h/2) and k r = n/Ar (n is the vertical mode number) 
since the radial and axial characteristic lengths are Ar and h/2, respectively. 

From Fig. [TTJ the Rayleigh's centrifugal instability, which occurs for all ( < when 
Va = and v = 0, is found to be suppressed by a strong magnetic field. This is consistent 
with Chandrasekhar 5^]. It is very interesting to see some growing modes when ( > 0, 
which is the MRI associated with the Stewartson layer. This instability also disappears 
with a sufficiently strong magnetic field. Comparing panel (a) with panel (b), we find that 
the magnetic field suppresses the modes with shorter wavelengths more strongly. This could 
explain why the Stewartson layer extends deeper into the bulk with a stronger magnetic field 
(Fig. [H]), by suppressing the growing modes with short wavelengths that would otherwise 
tend to smooth the velocity gradient. 

Viscosity also has a stabilizing influence (Fig. [l~2l . Given A = 1.5, the Stewartson layer 
is found to be stable if Re < 600 regardless of the vertical mode number. This could explain 
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Figure 12: MRI linear growth rates (7 s _1 ) as a function of Reynolds number Re and vorticity 
parameter £ at the junction between rings, = (ri+r2)/2 with A = 1.5. Parameters as in Tabled! 
Panel (a) the vertical mode number n = 0.25; panel (b) n = 1. 



why the Stewartson layer penetrates deeper into the bulk with increasing Reynolds number 
if the layer is steady (Fig. [10] (b), Re < 400). However at even larger Reynolds number 
the instabilities would destabilize the layer and presumably smooth it out except near the 
endcaps (Fig. [10] (b), Re > 400). If we could perform simulations without magnetic fields at 
the experimental Reynolds number (Re > 10 7 ), we expect that the profile of the azimuthal 
velocity v.s. radius at ~ 1 cm above the bottom endcap would be a tiny hump with a large 
slope of the azimuthal velocity \dv v /dr\ so that our experimental measurement would not 
resolve it, which matches the experimental evidence that an obvious Stewartson layer is not 
observed in a purely hydrodynamic experiment 21] . Unfortunately the current ZEUS code 
cannot afford a simulation with Reynolds number as high as the one in the experiment. 
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